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ABSTRACT 


A basin-scale acoustic tomography experiment was conducted in 
the northeast Pacific from May 1987 to September 1987. In this 
thesis, the stability of the forward model is analyzed. There 
are large non-linearities in the changes in travel time 
between ray paths for the four seasons. I constructed a model 
in which the change in warming in the upper 100 m of the ocean 
was due only to changes in surface solar irradiance. The value 
of the surface solar irradiance anomalies necessary to cause 
the tomography results for warming (Spiesberger and Metzger, 
1991) was computed. This value was larger than the actual 
value of surface solar irradiance anomaly which was computed 
using inputs measured by satellite (Chertock, 1989). 
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Chapter I 


Introduction 
1.1 Motivation and Overview 


Measurement of basin scale changes in ocean temperature, at 
the surface and in the interior, iS an important step in 
detecting global climate changes. Long range ocean acoustic 
tomography is an important tool for measuring the internal 
variability (Spiesberger and Metzger, 1991; Spiesberger et 
al., 1992). Satellites are an important tool in measuring 
basin-wide variations of temperature and radiation at the sea 
surface (Stewart, 1985). This thesis compares changes in 
basin-wide temperature in the interior with changes at the 
surface. I compute the surface solar irradiance anomalies 
necessary to achieve the topographically derived warming rate 
anomalies in the upper 100 m of the ocean along the experi- 
ment's path. The computation is based on a model in which the 
change in warming is due only to changes in surface solar 
irradiance anomalies. This thesis also investigates the 
stability of the models ee trace sound through the ocean. 

This thesis provides an introduction to tomography and to 
surface solar irradiance at the oceans surface in chapter I. 
Chapter II describes the forward model, examines the stability 
of the model from season to season, and compares the model to 


the experimental data. Chapter III examines the estimate of 





the surface solar irradiance anomalies and how it is computed 


Results and conclusions are given in chapter IV. 


1.2 Acoustic Tomography 


It is difficult to obtain time and space data coverage of 
the ocean by most traditional observational methods. Shipboard 
sensors can only measure one point of an observational area at 
a given time, and must move around the area taking many 
measurements to get good spatial data coverage. This can 
provide poor temporal resolution because the parameters being 
measured change during the ship's movements. Point measure- 
ments made by moored sensors provide good temporal resolution 
of one point over time, but providing sufficient areal 
coverage requires many such moorings. 

Ocean acoustic tomography was first proposed by Walter Munk 
and Carl Wunsch(1979) for measuring mesoscale(~100 km) proces- 
ses. It is a method of determining oceanic structure and 
variability through the inversion of acoustic travel-times to 
determine perturbations in sound-speed and currents and thus 
perturbations in temperature and density. Ocean acoustic 
tomography experiments are conducted by placing acoustic 
sources and receivers in the ocean. Each source transmits a 
pulse which is received at each receiver. The travel-time for 
each ray path is compared to the computed value for a model of 


known sound-speeds and the system is then “inverted" to 





estimate the sound-speed anomalies. 

One of the big advantages of tomography is that relatively 
few instruments are required to get good areal coverage of a 
large ocean area. If S sources and R receivers are placed 
around a volume and if P multipaths can be identified between 
them then there are S xX R x P pieces of data per measurement. 
If point measurements are used then there is only one datum 
per instrument(S + R data) at one time. Because sound travels 
at about 1500 m s’'! in the ocean the time to collect the 
tomography data is relatively short and thus can give a three 
dimensional synoptic measurement of oceanic parameters. 

Ocean acoustic tomography was successfully used ina 1981 
acoustic tomography experiment(Cornuelle et al., 1982; 1985) 
to map sound-speed anomalies on a 300 km square area near 
Bermuda. A 1981 experiment (Spiesberger et al., 1983) measured 
Gulf Stream meanders. The 1987 experiment by Spiesberger and 
Metzger (1991) and the 1983-1984 experiment by Spiesberger et 
al. (1992) demonstrated the use of sound to observe tempera- 


ture fluctuations over ocean basins (to 4000km). 
1.3 Surface Solar Irradiance 


The short wave radiation from the sun that is absorbed by 
the ocean, surface solar irradiance (Q,), 1S an important part 
of the heat balance of the atmosphere and the ocean. The other 


terms in the balance are latent and sensible heat fluxes and 





outgoing long-wave radiation. The long-wave radiation from the 


ocean, Q,, can be estimated using the formula (Gill, 1992), 


QO, = 0.985073 (0.39 - 0.05e,/) (1 - 0.6n2) (ae 


Where G 1S Stephen's constant for a black body, T,(°K) is 
temperature, 2 ust ells vapor pressure of water at the standard 
height(mb), iat is the fraction of sky covered by clouds, and 
the last factor,(1 - 0.6n¢), is a correction for clouds. The 
total upward flux of heat, Q, from the ocean is the sum of the 


fluxes of the individual processes, namely (Gill, 1982): ° 


Osage, * LE + Oe: (1.2) 


where Q, is the surface solar irradiance, for our purposes, 
computed by Chertock's algorithm(1989). Q. is the upward 
sensible heat flux, and LE the upward latent heat flux, with 
E the evaporation rate, and L, the latent heat of vaporization 


of water given by 


fi S 2 Ae yp eer Gis.) 


A simple model is used in this thesis to describe warming 
anomalies in the upper 100 m along the experiment's path. The 
model considers surface solar irradiance anomalies as the only 
source of warming anomalies. The other terms in the heat 
budget are not considered, including advection. I will use 
this model to quantify the magnitude of the surface solar 
irradiance anomaly necessary to cause the warming observed by 
tomography. 
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Chapter II 


Travel-Time Changes 
2.1 The Experiment 


The experiment used three acoustic sources and seven receiv- 
ers. The sources were deployed in a triangle, approximately 
1000 km on a Side, north of Hawaii (Figure 2.1). Source (S1) 
was deployed in about 5500 m of water on a taunt mooring at 
667 m depth. It transmitted every 2 hours, 12 times a day, 
every fourth day. The pulse-like signal's acoustic frequency 
was centered at 250 Hz with an rms bandwidth of 63 Hz. Thus it 
had a pulse resolution of (63Hz)"' * 16 ms (Spiesberger and 
Metzger, 1991). The source level was 191 dB re 1 pPa @lm. The 
time of the source transmission of the signal is known to 
about 1 ms with the use of a rubidium frequency standard. The 
time was known to the same accuracy at the receiver. The 
geographic positions of the sourceS were measured to an 
accuracy of about 10 m with the Global Positioning System. 
Mooring motion was tracked with the use of four bottom-mounted 
transponders and an A nptureveratareieeare on the mooring near the 
source (Liberatore, 1985). The receivers are bottom mounted 
off the west coast of the United States. Their positions are 
known to an accuracy of about 120 m (Spiesberger and Metzger, 
1991). 


The Acoustic path I studied extends over about 2000 km and 
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was insonified for 120 days from May 1987 to September 1987. 
It was necessary to raise the SNR (Signal to Noise Ratio) 
because the source could not transmit a pulse powerful enough 
to be heard above the background noise level. The SNR was 
increased about 36 dB when the received acoustic Signal was 
complex demodulated, low-passed filtered, correlated with a 
replica (phase-only-filtering), and averaged over four 
sequence periods (Spiesberger et al., 1989). The data were 
corrected for mooring motion and clock drift to within the 
limits noted above. The source used for this experiment is 
designated S1 and the receiver is designated Rl on Figure 2.1. 
Source S3 and receiver R2 were used by Spiesberger and Metzger 


(1991). 


eZ 





© Acoustic Source 

CL) Acoustic Receiver 
Hydrographic Stations: 

A 1Aug- 14Sept 

4 14 Sept - 24Sept 





240 °E 
Figure 2.1 Schematic diagram of the experiment with the posi- 
tion of the sources and receivers indicated. The acoustic path 
used in this paper is indicated with the heavy dashed line. 
From Spiesberger and Metzger, 1991. 
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2.2 Ray Tracing 


The speed of sound in the ocean depends on temperature, 
pressure, and salinity. Del Grosso's (1974) algorithm for 


sound-speed in sea water is, 


Cee ae tee, + Ne. + NC, + ACen, (2.1) 


where, 
Coop = 1402.392 


mere 5 Ono — 550946410717? + .2215153X10-T? 
Wee = 1.329525 + .128955X10--S* 
Bien = .156059P + .244998X104*P* - .883392X10-'P? 


opr sie x Or ro + ,635X10°TP © .2654X10 'T*P* 
—-.159xX10°TP* + .522X109TP?-.438X10-6T?P 


smn somo nt 960.1 Osr > + .485X10 -TS*P 
-,.340X10°°TSP 
where T is temperature in degrees Celsius, S is salinity in 
parts per thousand, and P is pressure in kilograms per square 
centimeter. Temperature and pressure have a much greater effect 
than salinity on sound-speed in most regions of the ocean. 
We will now develop a solution to the acoustic wave 
equation in terms of ray theory. The wave equation is (Tolstoy 
and Clay, 1966) 


V2 -_ ah Gag =o) Zee 
zs Camale - ) 





where V is the Laplacian operator, p is the acoustic pressure 
and c is the sound-speed, a function of salinity, pressure, and 
temperature which may vary with spatial coordinates(i.e., ¢c = 
c(x,y,2))-. For a harmonic source, with w being the angular 
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frequency, ea tHt the wave equation becomes the Helmholtz 


equation: 


Wine Kp = 0 (223) 


with k being the wavenumber. Rewriting the wavenumber k as 


Cc 
ee ee (ania) 
0 


ene. (2.5) 


is the index of refraction. I write acoustic pressure in the 


form 


p= Ae tkoS (x.y. z) (2.6) 


where A is the wave amplitude, a function of (x,y,z), and k,& 
is the phase of the wave. Substituting equation (2.6) into 
equation (2.3), and collecting together the real and imaginary 


parts results in the equations 


VA - kfA[(n? - (VS)2] =0 (2.7) 


and 


2VA:' VS + AV’S = 0 (2.8) 


If we make the assumption: 


VA 
k*A 


<i (2.9) 











1-5 





then this implies the rate of variation of the wave phase, per 
wavelength, of the vertical component of wavelength is small 
(Tolstoy and Clay, 1966). Equation (2.10) can be shown to 


follow from equation (2.9) (Frisk, 1991) 


ac ey 
pei ad 2.10 
| < A 


where A is the wavelength. This is a necessary, but not 
sufficient, condition for the ray acoustics approximation. It 
is also necessary that the reference point is not in the shadow 
zone or near its boundaries, and also not near caustics 
(Brekhovskikh and Lysanov, 1982). Even simpler, the medium must 
vary slowly over an acoustic wavelength. This is the geometric 
approximation. 


Applying the condition (2.9) to equation (2.7), results in: 
(VS\e = ne (2.11) 


or, expanding the gradient operator: 
(52) ‘ (se) + (52) ae (2.12) 

Ox Oy OZ 
This is the eikonal equation, which forms the basis for ray 
theory. A physical picture of rays follows directly from the 
eikonal Beer ion. Surfaces of constant phase (wave fronts) are 
given by S = constant, and the lines orthogonal to the wave 
fronts (VS) define the geometry of the rays(Brekhovskikh and 
Lysanov, 1982). The rays represent the paths along which 
acoustic energy propagates. The amplitude of the rays is given 
by equation (2.8), the transport equation. For our purposes, we 
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consider propagation in the vertical plane and let sound-speed 
be a function of depth only (1.e., c = c(z)) between sound- 
speed profile segments; a good approximation because of the 
vertically stratified nature of the ocean. Thus, if an acoustic 
signal has a wavelength that is much less than the length scale 
being examined then the geometric approximation is valid and 
ray theory can be used to model the acoustic field. 

In middle latitudes the sound-speed decreases from the 
surface to a minimum axis at ~1000m due to temperature decreas- 
ing, and then increases to the bottom due to pressure increas- 


ing in the nearly isothermal depths of the ocean(Fig 2.2). 


“ 
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= 
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—— 
Q 
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000 
1475 1480 1485 1490 1495 1500 1505 1510 1515 1520 1525 


Sound-speed in m/sec 





Figure 2.2 A typical sound-speed profile. 
Sound-speed is computed with Del Grosso's 
alg. using Levitus data (Levitus, 1982). 
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Sound rays bend towards minimum sound speed according to 


Snell's law of refraction 


Gosu) cost) 


Cy Co 


= Constant (2.23) 





where c, is sound-speed in the layer and 6, is the angle the ray 


makes with the horizontal. 


Sound 
Speed 





Acoustic 
Ray 


Figure 2.3 Snell's law of refraction for c, < c,. 


This shows that sound rays will constantly refract toward the 
area of lower sound-speed. Because there is a minimum in sound- 
speed, at about 1 km depth at mid-equatorial latitudes, sound 
energy is bent away from the high loss regions of the surface 
and bottom. Each ray path samples a different vertical section 
of the ocean depending on the launch(departure) angle from the 
source. Travel times for ray paths are different because sound 


speeds are different along their paths. These times are given 


by 
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qT = {_ds_ (2.14) 
/ c(s) 


where T is the travel time along the ray path for a particular 
ray, p is the ray path, ds is the arc length, and c(s) is the 
sound-speed at each point along the path. The ray paths are 
such that only a few rays actually connect the source and the 


receiver; these are defined as eigen-rays. 


2.3 The Ray Tracing Program 


In this project, we used the range dependent ray tracing 
program MPP (Multiple Profile ray tracing Program) developed by 
C. W. Spofford. The sound-speed field is linearly interpolated 
in both depth and range in specified triangular sectors. The 
sound-speed field is continuous everywhere, but is's gradient 
is discontinuous at the triangular boundaries. The bottom 
bathymetry is represented by piece-wise linear segments. Output 
of the program includes eigen-ray arrival times and transmis- 
Sion loss (calculated from geometrical spreading and losses due 
to boundary reflections), along with a history of the eigen-ray 
trajectories. 

The geodesic path between the source and receiver was 
computed using the WGS84 reference ellipsoid of the Earth 
(Defense Mapping Agency, 1987). We compute sound speed profiles 
from the Levitus seasonal data base(Levitus, 1982) using Del 
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Grosso's algorithm (equation. 2.1). Bathymetry points along the 
geodesic are from a digital data base (National Geophysical 
Data Center, 1987). Using a flat-Earth coordinate transforma- 
tion (Ben-Menahem and Sarva, 1981), rays are traced in a 
Cartesian coordinate system. 

MPP computes two types of rays (Spiesberger et al., 1991). 
The first type, the geometric ray, passes through both the 
source and receiver. The second type, the diffracted arrival 
(Brown, 1982; Pierce, 1981), does not pass through the receiv- 
er, but through points near the receiver. Because propagation 
is at finite frequency, energy leaks into the Beever by 


@iffraction as an exponentially decreasing wave. 
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2.4 lLinearizing The Forward Problem 


The travel-time of the i1,, ray along path [l, is given by 


ds 
T, = <n 2. 
: 5 c(s) + uct C2on 2) 


Peme sds 1S 4 differential arc-length along the path T., c(s) 
is the sound-speed, u is the current component along the ray, 
° is a dot product, and 7 is a unit vector tangent to the ray. 
The travel time of a given ray 1s dependent upon the path 
length, sound-speed, and current velocity along the ray path. 
Variations in sound-speed and current will lead to Pee tilions 
in travel-time and a change in the ray path. Hamilton et al., 
(1980) showed that there is a small change in travel-time 
associated with this change in path length if sound speed and 
current variations are small. However, the change in path 
length leads to a different sampling of the oceanic mediun. 
Now we will examine the relative size of the terms in the 
denominator of the integrand of equation (2.15). A typical 
value for current speed is u = 0.1 m/s and a typical sound- 


speed is c = 1500 m/s, so 





4 =0(1074) <1. 
‘alll 


Typical values for the vertical shear of current and sound 





du Oem s i dc 5m/s 
d — = —————_. = 0(107 — = =O (110 2 
spee are —— IGG ( ) and cP Ware ( ae 
ac. du 
SO FE 1s usually one order of magnitude larger than rea 


The refraction of rays is dominated by the sound-speed gradi- 


ent, and the current can be ilgnored in ray tracing, to first 


PAA 





Order. 
The path of the ray in equation(2.15) is a nonlinear 
function of the sound-speed field, so we linearize about a 


reference state, c,, as in: 
Gy = "Cny) + OEY) (2.16) 
where 6c is the perturbation from the reference field. Also 


[5c(x)| < c,(x) (2.17) 


Substituting (2.16) into the travel-time integral (2.15) gives 


ale i CS a 
De i ERGO Cc (S) at BoD) 


Expanding the integrand, keeping only the leading order terms, 


and dropping the current velocity term yields: 


Tes f oe OCS) aS (2.19) 
ae eres) r,; c2 (S) ° 
Thus the travel-time perturbation can be expressed as 
6c(S) ds 
eT, = - ( 2e(S) ds | (2.20) 
Lh Tey 


Equation (2.20) is linear in the sound-speed perturbation, éc, 
if [, is assumed not to change as the sound-speed field 


perturbs. The forward problem is finding the 6T, when éc is 


Known. 


Sound-speed is directly proportional to temperature (6), 
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with an empirical relationship given by: (Del Grosso, 1974) 


6C _ 459 (2.21) 


Co 
where a = 3.2 X 10%. Thus equation (2.20) can be considered a 


linear relationship between the travel-time perturbations and 


perturbations in temperature. 
2.5 Identification of Eigen-rays 


To identify eigen-rays for the spring and summer seasons I 
computed a daily average of the 12 arrival records for eh aan 
data was available as prescribed by Spiesberger et al., (1980). 
I compared the output of MPP for spring and summer seasons to 
the daily averages. The output of the MPP ray tracing program 
is shown in Figures 2.4a and 2.5a. The Levitus data base for 
the spring and summer seasons is used to generate sound speed 
profiles with Del Grosso's algorithm(equation. 2.1). The 
Levitus data base is divided into the four seasons according to 


Table 2.1. 
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Table 2.1 Levitus data base seasons 


Calendar Months Calendar Days 


‘| = 4 
Spring May - June 121 - 212 


August - october | 213 - 304 
November - January | __305 - 032 
February - April | 032 - 120 












Daily averages for Days 189 and 257 are shown in Figures 
2.4b and 2.5b. Six eigen-rays (all geometric ray arrivals) ane 
identified for the spring season. The last pair of arrivals, in 
Figure 2.4b, are predicted as diffracted arrivals from MPP 
output. Because there are a large number of diffracted rays 
from MPP output grouped together near the peak arrival time, 
the individual diffracted rays could not be resolved with 16 ms 
resolution. The individual rays of the groups near the last two 
arrivals have different paths, sampling different parts of the 
ocean and thus are not used here for tomographic inversions. 

For the summer season seven geometric eigen-rays were 
identified. The last arrival was predicted from MPP as a 
diffracted arrival and the individual diffracted rays could not 
be resolved. The geometric rays that have larger travel times 
than the identified eigen-rays intersect the bottom at large 
angles and are not observed at the receiver. 


The daily average for day 189 is compared with the MPP spring 
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ray trace output in Figure 2.6a. Figure 2.6b shows the daily 
average for day 257 compared to the MPP spring ray trace 


output. 
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Figure 2.4 (a. Top) MPP output for the spring Levitus data 
base. (b. Bottom) Daily average for julian day 189. X-axis is 
time(s), Y axis is relative magnitude and the start time is 


about 20 minutes. 
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siv4 Forward geometrics & difracted rays from mpp, Summer only 
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Figure 2.5 (a. Top) MPP output for the summer Levitus data 
base. (b. Bottom) Daily average for julian day 257. The X-axis 
is time(s), the Y-axis is relative magnitude, and the start 


time is about 20 minutes. 
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Figure 
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m2 Daily average for s1v4 for day 257 
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2.6 (a. Top) Day 257 average overlaid with MPP summer 
(b. Bottom) Day 189 average and MPP spring output. The 
is time(s) and the Y_axis is relative magnitude. The MPP 


was shifted 60ms. 
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Table 2.2 shows the difference in travel time (ms) between 
the observed daily average peak arrival time for days 189 and 
257 and the predicted eigen-ray travel time for the spring and 
summer seasons respectfully. Predicted times were less than 
their corresponding observed arrival times for all rays (Table 
2.2). A time difference of 0.1 second between measured and 
predicted travel times in the upper 1000m along our ray paths 
equates to less than 0.1°C temperature difference (Spiesberger 


Baal., 1983). 


Table 2.2 Daily average peak arrival time (measured) - ray 
travel time (predicted); for days 189 and 257, Spring and 
Summer. 


Peak/Ray Spring travel time Summer travel time 
difference (ms) difference (ms) 
measured - predicted measured - predicted 





Not Applicable 


Zo 





2.6 Variations in the Forward Model 


I ran the MPP ray tracing program for the four seasons using 
the Levitus climatological data base of the oceans to compute 
sound speed profiles with Del Grosso's algorithm (equation. 
2.1). We identified six eigen-rays (all geometric arrivals) for 
the spring, fall, and winter seasons and seven for the summer 
season. The first 100 km of the propagation path of the six 
elgen-rays common to all four seasons are shown in Figure 


2.7(a-f). 
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eigen ray #2 path 






solid = spring, summer = dashed, fall = dotted, winter = stars 


eigen ray #4 path 
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Figure 2.7(a —- d) Paths of first 4 eigen rays for the lst 100 
km. X-axis is horizontal range(km) and the Y_ axis is depth(m). 
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S1V4 eigen ray #6 path 
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Figure 2.7(e - £) Paths of last 2 eigen rays for the lst 100 
km. The X-axis is horizontal range(km) and the Y-axis is 
depth (m) 


To compare the nonlinearity of the forward problem for the 
four seasons we follow ina similar fashion to a derivation 
found in Spiesberger (1985). The reference travel time between 


a fixed source and receiver along the reference ray path IT, is 


given by: 
ies ff 2 (2.22) 
ee oF 
where Cc, = c)(X,Z) 1s the initial sound speed. The perturbed 


travel time along the perturbed ray path ©, is given by: 
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= ds 
jae aca r 7 (2.23) 


where dc = 6c(x,z) is the perturbation in the sound speed. The 
perturbed ray path becomes the unperturbed ray path as 6c > 0. 
The change in travel time between the initial and perturbed 


state is given by: 


6T, = T, - T= eis aS 


2.24 
hee, hoc If, GC, ( 


Because |Sc/c,| * -003 for mid-ocean sound speed perturbations 


we can expand T, as: 








T, = w= f 8 - 26 « Oey a5. | Rees) 
= dik, Gy ayes paren e “on 
Thus 6T, becomes: 
ve 2 OC) ooo | 27 BE TARAS 
man veer fer on aeons 
or 
2 
ee | - feito - 22... . (2.27) 
eaten roerer re, cé C: 





This can be written as: 


6c 
Co 


a, = ( SB - (Bl - fas 3th - 


*F _ (2.28) 
Yr, Co re) Co Co 





the last term of equation (2.28) can be written as 
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1- S24---Je fasts + an (2.29) 
Co Np 





- fos 23 


The first term of on the right-hand-side of equation (2.29), 


6 
6T, = -| ds—= (2.30) 


is a linear approximation for the exact travel time change. The 
most important part of the linearizing is that ©, =[T, (i.e. the 
perturbed path is the same as the reference path). 


Defining 


Hilly ol Cp (2.31) 
Tr, Cg Try Co 


and writing 6T, from equation (2.29) as 


éT, = ass - fase + [ be _ (2.32) 
Sy] es 1” (eee MmGgic. + OC) 
allows the exact travel time change to be written as 
ore = OL, + OT, + OT, « (2.33) 
This can be expressed as 
6T, = T, - T, = 6T, + NL (2.34) 
or 
eae Le) = Ol, (2.35) 
where 
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Picola OL, (2.36) 


NL contains all of the nonlinear terms of the forward problem. 
If the initial ray and the perturbed ray coincide everywhere 
then NL = 0. The percentage of nonlinearity of the forward 


problem is: 





NL 
2 24530 
$ NL 109 | ( ) 


NL is a measure of the nonlinear effects, primarily due to 
changes in the ray path, it is a measure of the errors of the 
linear approximation 6T, verses the exact travel time Shands 
5T,. 

For our experiment we used the Levitus summer sound speed 
profiles and the MPP output ray paths as our reference state 
and the other seasons as our perturbed states. T, and T,, the 
travel times along the reference and perturbed paths are 
obtained from MPP output. *NL was computed using the summer as 
reference for all six of the eigen-rays for each of the other 
seasons. The results are summarized in Table 2.3. 

Most of the rays have large nonlinearity, thus 6T, is not a 
good approximation to the araee travel time changes for these 
rays. There are exceptions; for example eigen-ray # 3 for two 
seasons and eigen-ray # 4 for all seasons have small nonlinear- 


mey . 
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Table 2.3 Seasonal variations in the forward model. 6T, is 
given by Eq. 2.24, 6T, by Eq. 2.30, NL by Eq.2.35 and eNL is 
given by Eq. 2.37. The reference paths are computed from the 
summer time values. 


| season | srs) | sts) | wt(s) 


Eigen ray # 2 


es 


a a > ae pees eee cul 6 
winter 


Eigen ray # 2 


Set 
Eigen ray # 3 


Elgen ray # 4 


— a eT 


joe 


Fall eS ra 

| ass] e302] =e | 27-98 
Eigen ray # 6 

[winter | e276] eta] tz | ts 

raul | —en90| cosa | tz | 63 


Winter mae .144 50.4 


| ee 
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Table 2.4 shows the first upper turning depths for the six 
eigen-rays for each season. 


Table 2.4: Eigen ray upper turning depth for the four seasons. 


Rf 
[spring | o | o | o | o | a5 | 23.2 | 
| summer | 16.6 | 22.1 | 30.9 | 32.2 | 38.9 | 40.8 | 
fra | o [| o [| o [| o | o | o | 
=i Ss a Se a 








2.7 Observational Variations in Travel Time 

Figure 2.8a is a gray scale bit plot of the observed data 
with the darker shades corresponding to stronger amplitudes. 
Peaks (arrivals) were picked from the data using a variable 
size window scheme that determined the largest signal within 
the window. The picked peaks are shown in white in Figure 2.8b. 
The daily average multipath arrival (peaks) trends are shown in 
Figure 2.9. Variations in the travel time of the daily averages 
are larger for the later arrivals than the early arrivals. 
Multipath arrivals start arriving earlier around year day 180 
(late June) with the sharpest decrease in arrival times between 
late June and late July. The later arrivals correspond to rays 
with smaller launch angles, thus sampling less of the deep and 


surface water. 
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Figure 2.9 The daily average multipath peak trends. The start 
times of the trends are offset by 0.2 s and initialize to zero 
at day 140. The Y-axis is in seconds and the X-axis is 1987 


yeardays. 
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Chapter III 


Surface Solar Irradiance Variations 


3.1 Satellite Data 


The four components of the surface heat flux budget are: 1) 
the net surface short wave radiation; 2) the net surface long 
wave radiation; 3) sensible heat; and 4) latent heat. Chertock 
(1989) generated a global seven year record of net solar 
irradiance at the ocean surface for the period between November 
1978 and October 1985. Her algorithm is based ona Se aa aeive 
transfer theory. Surface solar irradiance is the net radiation 
from the sun entering the surface of the ocean. The surface 
solar irradiance was computed by subtracting the sum of solar 
energy absorbed by the atmosphere plus the solar energy 
reflected by the ocean-atmosphere system from the solar energy 
at the top of the atmosphere. 

The solar energy absorbed by the atmosphere is a modeled 
quantity computed by Chertock (1989). Albedo is the ratio of 
the radiation reflected by a body to the amount incident upon 
it. Planetary albedo is measured by satellite for this experi- 
ment. The solar energy at the top of the atmosphere is a known 
quantity calculated according to the astronomical sun-earth 
geometry. 

The planetary albedo used as input for Chertock's algorithm 


came from the Nimbus-7 satellite. The Earth Radiation Budget 
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(ERB) experiments Wide Field Of View (WFOV) radiometer on 
Nimbus-7 provided the measurements (NASA, 1984b). Twenty-two 
spectral channels of radiation were measured by the ERB 
experiment instruments. The infrared radiation emitted by the 
earth and the solar radiation reflected by the earth are 
measured by the four fixed sensors of the WFOV radiometer part 
of the ERB package. The footprint (ground resolution) of the 
WFOV instrument aboard Nimbus-7 is approximately 1000km (NASA, 
1984a). The WFOV radiometer measures the reflected solar 
radiation in the spectral range between 0.2 and 3.8um (NASA, 
1984b). | 

The Nimbus-7 satellite was launched on October 24, 1978 into 
a nearly circular, sun-synchronous orbit at an altitude of 
950km. The sun-synchronous orbit provides two passes per day, 
near local noon and midnight for equatorial crossings. Measure- 
ments by a satellite in a sun-synchronous orbit are always made 
at the same times for a given location (Stewart, 1985). This 
produces a biased distribution of data over the diurnal cycle. 
If the parameter being measured varies diurnally, as planetary 
albedo does, the bias must be accounted for in the computa- 
tions. 

The diurnal variations in the albedo are quite large, a 
Bacco: Of 1:5 to 2 for clouds and a factor of 4 for ocean 
surfaces (Hucek et al., 1987). There are also diurnal varia- 
tions in cloud cover (Minnis and Harrison, 1984) that add to 


the diurnal variability of the albedo. 
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Albert Arking of the Goddard Space Flight Center designed a 
global reflectance model to convert the measured WFOV planetary 
albedo to a daily average planetary albedo. He used the Nimbus- 
3 directional reflectance models designed by Raschke et 
al.(1973) to estimate that 90% of the solar radiation reflected 
came from cloud-land reflectance and about 10% was from ocean 
reflectance (Kyle et al., 1986). Table 3.1 shows the Nimbus-3 
directional values derived by Raschke et al., (1973) and the 


Nimbus-7 ERB WFOV composite-directional reflectance model. 


Table 3.1 Normalized directional reflectance used in Nimbus-7 
ERB algorithms to calculate planetary albedo diurnal correction 
factors (Kyle et al., 1986) 


| zenith angle land Composite* 






1.56 
1.78 
1.60 1.34 | 


* Composite is 90% cloud/land plus 10% ocean. 






The planetary albedo data was processed using calibration 
adjustment procedures detailed in Kyle et al.(1985). Daily mean 


planetary albedo data is computed (Kyle et al., 1986) from the 
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instantaneous measurements of the Nimbus-7 ERB experiment. 
To use the computed daily average planetary albedo for 
studies of the global radiation budget it is necessary to 


weight the measured solar flux by the daily average insolation 


at the top of the atmosphere (Hucek et al., 1987). A,,,, the 
monthly mean planetary albedo, is computed from 
M => 
YK F(o,4, €3) 
A fig (3.1) 


sat 


PreeslOy Moe ep 


Ms 


a) 
WW 
be 


Where M is the number of days in the month and k, is the albedo 
normalization constant for the i,, day. F is the instantaneous 
femrbratced measurement of the reflected solar radiation. F,,, is 
the amount of solar radiation that would have been reflected by 
the incident object if it had a perfectly reflecting diffuse 
surface. F and ee adncwene saviekagc .OLrar sand bi For all the 
measurements within the field of view for a single pass of the 
satellite. ¢ is the latitude of the sub-point target area, 4A is 
eiemrongitude of the sub-point target area, and t, is the day 


of the year. 
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3.2 Computation of Surface Solar Irradiance 


The solar irradiance at the ocean surface was computed by 
Chertock using her algorithm based on radiative transfer theory 
(Chertock, 1989; Frouin and Chertock, 1992; Chertock et al., 
1992). Chertock developed a layered model of the atmosphere 
based on plane-parallel theory to compute surface solar 
irradiance. The top layer is a clear sky atmosphere which is 
above, and decoupled from, an effective cloud layer. The solar 
irradiance at the surface is a result of the radiation passing 
through the clear sky and cloud layers. The solar Santon 
incident at the top of the atmosphere is reduced by scattering 
and absorption through the clear sky layer by aerosols and 
molecules. The solar radiation is further reduced as it passes 
through the cloud layer by absorption and multiple reflections. 
Chertock used this algorithm (Chertock, 1989; Frouin and 
Chertock, 1992) to convert the monthly mean planetary albedo 
data provided by the satellite to surface solar irradiance. 

The largest contributors to atmospheric absorption of 
radiation in the solar spectrum are oxygen (0,), carbon dioxide 
(CO,), Ozone (03), and eee (H,0) (Wallace and Hobbs, 
#977). The first two, O, and CO,, are well mixed in the atmo- 
sphere and can be considered to be at constant concentrations 
for radiation study purposes. However, the concentration and 
distribution of ozone and water vapor are highly variable, both 


in time and space, thus having a variable effect on the amount 
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of reflected radiation measured by satellite (Wallace and 
Hobbs, 1977). 

Chertock used the algorithm created by Tanré et al., (1986) 
to compute atmospheric absorption and scattering in the clear 
sky layer. This code, referred to as "5S" (Simulation of the 
Satellite Signal in the Solar Spectrum), relates the radiative 
properties of the planetary system to the satellite measured 
solar radiance in clear atmosphere. The "5S" code was used to 
compute the diffuse atmospheric transmittance due to scattering 
by aerosols and molecules. The "5S" code also calculated the 
transmittance due to gaseous absorption. The code also Sores 
the portion of the measured flux that is from radiation 
reflected by the atmosphere and returned to space without ever 
reaching the cloud layer. The 5S code uses the following 
inputs: 1) the total water vapor and ozone amounts, and the 
vertical distributions of the ozone and water vapor; 2) the 
geometrical conditions; 3) The ground reflectance (a function 
of wavelength); and 4) an aerosol model specifying maritime or 
continental background and aerosol concentration (a function of 


wesibility) . 
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3.3 Anomalies of Surface Solar Irradiance 


Chertock (1989) produced a seven year (November 1978 - 
October 1985) record of net solar irradiance at the ocean's 
surface using a numerical code to implement her algorithm. I 
provided her with Nimbus-7 ERB WFOV planetary albedo data 
(obtained from NASA's Climate Data System User Support Office 
at NASA Goddard Space Flight Center) covering the months of 
interest for the Pacific experiment. She used these data to 
compute monthly mean surface solar irradiance values for April 
1987 through September 1987. The surface solar irradiance ane 
are on a 9° latitude-longitude spacial grid box over the 
world's oceans. 

Climatological monthly means were computed by averaging the 
data for the reference period (November 1978 - October 1985) 
excluding El Nino years (January 1982 - December 1983). Errors 
in the surface solar irradiance values are estimated to range 
between 10 and 20 W/m* (Chertock et al., 1992). Most of the 
anomalies shown in Figure 3.1(a-f) are within the error limits 
and thus are not significantly different from zero. The anoma- 
lies of surface solar irradiance were computed by subtracting 
the climatology value from the monthly value. The variation 
from the monthly means for the months from April 1987 to 
September 1987 were computed by Chertock for each grid box in 
the Pacific experiment area. The values shown on Figures 3.1(a 


- f) represent averages for the entire nine degree grid box. 
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APRIL 1987 SURFACE SOLAR IRRADIANCE ANOMALIES (W mes) 





166.5°W 157.5° 148.5° 13 De * T3055 Ls 


Figure 3.1a April 1987. The X-axis is longitude and the Y- 
axis is latitude. 
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Figure 3.1b May 1987, Acoustic path between S1 and Rl is 
indicated with a dashed line 
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JUNE 1987 SURFACE SOLAR IRRADIANCE ANOMALIES (W m’°) 
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Figure 3.1ic June 1987, The acoustic path analyzed by Spies- 
berger and Metzger (1991) S3 - R2 is indicated with a solid 
line. 
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Figure 3.1d July 1987 
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AUGUST 1987 SURFACE SOLAR IRRADIANCE ANOMALIES (W m’‘°) 
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Figure 3.1le August 1987 
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Figure 3.1f September 1987 
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To examine the possible effects of changes in surface solar 
irradiance on acoustic travel times I will look at heating 
rates in the upper 100 m as Spiesberger and Metzger (1991) did. 
I will use a model in which the only parameter affecting the 
heating rate in the upper 100 m is changes in surface solar 
irradiance. The model does not consider changes in other heat 
budget terms (section 1.3) nor does it consider advection. 

To look at changes in warming, I first determined the 
average surface solar irradiance anomaly along the acoustic 
path during each month of the experiment (the dashed line in 
Figure 3.1(b)). I also determined it for the acoustic path 
analyzed by Spiesberger and Metzger (1991) (the solid line in 
Figure 2.1 and in Figure 3.1(c)). The averages were determined 
by weighing the value in each grid box by the percentage of the 
path that went through the grid box. The average monthly 
surface solar irradiance anomalies along both paths are listed 
in Table 3.2. Note that the magnitude of the anomalies are 
primarily within the error limits of the surface solar irradi- 


ance data noted previously (10 W/m? to 20 W/m‘). 


Table 3.2 The average surface solar irradiance anomalies (w/ 
m*) along the acoustic paths for the sections S1-Rl and S3-R2. 


ay | gune | uty | august | sept. 





I computed the change in average monthly surface solar 


DAL 





irradiance anomaly from one month to the next, 6A. For example, 
for August to September for S3 - R2 path, 6A = -17.7 - (- 


ae w/in = 6.6 W/m* (Table 3.2). The values of 6A are between 


1.6 and 7.2 W/m*. The standard deviation of 6A is /2 * standard 


deviation of each monthly anomaly; namely 2 ¥*20W/m?=29 . Thus 
values of the differences of the surface solar irradiance 


anomalies are not significantly different than zero. Based on 
this model, there is no significant change in the warming rate 
from the historical mean if surface solar irradiance anomalies 
is the only cause of changes in warming. 

Spiesberger and Metzger (1991) used tomography to show that 
the warming between August 14 and September 14, 1987 was less 
than normal by 0.2°C+t0.056°C in the upper 100 m along the path 
from S3 to R2 (Figure 2.1). If I assume that this change in 
warming was due only to changes in surface solar irradiance, I 
can compute the surface solar irradiance anomaly decrease over 
the month necessary to achieve this decrease in warming. 
Although not realistic, this assumption allows me to determine 
the magnitude of the surface solar irradiance anomalies 
necessary to change the warming by the above amount. Using the 
formula from Wyrtki and Uhrich (1982) to compute 6T, the change 


in warming(°C) is: 


Or, 
56T = 5 CD (3.2) 





where p is the density of water (10° gm/n’), Cc, is the specific 


OZ 





heat capacity at constant pressure (4.2 J/gm °C), D is depth 
(100 m), and 6H is the change in heating (J/m*). For 6T = 

-0.2°C + 0.056°C this would result in 5H = - 8.4 x 10% J/m* + 
2.35 x 10’ W/m*. To get this change in 6H, I computed the 
necessary change in surface solar irradiance anomaly (6A) for 


30 days with, 


_ 6H(J/m*) . 1lday .1month 
6A = Oe (3-3) 
3600s/hour 24hour 30day 


Equation 3.3 yields -32.4 W/m* + 9.1 W/m* for a month. To get 
the results of Spiesberger and Metzger using only changes in 
the monthly surface solar irradiance would require Se ieeerence 
of -32.4 W/m? in the monthly anomalies from one month to the 
next. 

For the above results I will consider two standard devia- 
tions away from each mean. This gives 6A from tomography of 

-32.4 W/m* + 18.2 W/m* and from the surface solar irradiance 
model 6A is 6.6 W/m* + 58 W/m*. These values are less than two 
standard deviations apart. 

I wish to restate that acoustic tomography measures a 
vertical slice of the ocean integrated along the acoustic path. 
Surface solar irradiance anomalies are a monthly average over 
a 9° grid of the ocean surface. The model considers only 
changes in surface solar irradiance and neglects changes in the 
advection and other heat budget terms. Changes in the neglected 


terms may dominate changes in warming. 
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Chapter IV 


Discussion and conclusions 
4.1 The Forward Model 


We used the Levitus climatological data base with the Del 
Grosso sound speed algorithm to generate sound speed profiles 
as input to our forward model. The MPP ray trace program was 
used as the forward model. MPP produces both geometric and 
diffracted ray arrivals as output. We identified six eigen-rays 
for the spring, fall, and winter seasons and seven RoERthe 
summer season. 

We compared the MPP output for the spring and summer to the 
1987 experiment data. The MPP forward model did a good job of 
predicting arrival times. The sound rays take about twenty 
minutes to travel between the source and receiver. MPP predict- 
ed their arrival time within 0.1 seconds(Table 2.2). Tempera- 


ture fluctuations can account for a 0.1 s offset. 
4.2 Variations in the Forward Model 


The seasonal stability in the ray paths of the forward model 
was examined by computing the non-linearity in the travel time 
changes as a function of season. It was shown (Table 2.3) that 
there is a large percentage of non-linearity in the changes of 


travel time (>30%) between the seasons for most of the eigen- 
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rays using summer ray paths as the reference; the summer eigen- 
ray paths are not a good approximation for the other seasons 


for tomographic inversions. 


4.3 The Observed Acoustic Data 


The travel time decreases between year-day 140 and year-day 
260. The change in travel time over the 120 day period is 
largest for the later arrivals(Figure 2.9). The later arrivals 
show decreases in total travel time which indicate faster ound 
speeds along the ray's trajectory. The later arrivals corre- 
spond to rays with smaller launch angles, that sample less of 
the near surface waters. Multipaths closest to the surface show 
much smaller changes in travel time. The earlier arrivals show 
an increase in travel time from year-day 140 to year-day 195 
and then a decrease in travel time until the end of the 
experiment, year-day 260. An increase in travel time corre- 
sponds to slower sound speeds along the ray path, indicating 


cooler water. 
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4.4 Surface Solar Irradiance 

The solar irradiance anomalies at the oceans surface, 
computed by Chertock, along the experiment's sound path range 
from -4 W/m* + 20 W/m* to -35 W/m? + 20 W/m*. The magnitude of 
the surface solar irradiance anomalies shown in Figures 3.1 (a- 
f) are typically within the error limits of the surface solar 
irradiance data and so the changes in warming due to from them 
have no statistical significance. 

I considered a simple model in which ocean temperature 
changes in the first 100 m of the water column are a function 
only of changes in the surface solar irradiance. The model does 
not consider changes in the other terms of the heat budget 
(section 1.3) nor does it consider advection. I computed the 
surface solar irradiance anomalies (using the model) necessary 
for the change in warming found by Spiesberger and Metzger. 
They derived the change in warming from tomographic inversions 
for the path they analyzed in the experiment area. They found 


below normal heating from August to September. 
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4.5 Conclusions 


Surface solar irradiance is only one of four terms in the 
heat budget, the other three being long-wave radiation, 
sensible heat, and latent heat. We must also consider advection 
of heat when looking at the ocean's heat budget. No direct 
relationship between surface solar irradiance anomalies and 
travel time changes can be inferred without doing a complete 
heat budget including advection effects. 

Acoustic tomography has been shown to be capable of measur- 
ing the thermal properties of the ocean on a basin eae 
including inter-annual temperature variations (Spiesberger and 
Metzger, 1991; Spiesberger et al., 1992). A long term tomo- 
graphic experiment covering the entire ocean, along with a heat 
budget computation including the terms discussed in section 1.3 
and advection, would be important in measuring changes in 


warming of the ocean. 
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